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I. PRELIM IHARY OBSERVATIONS 



The disturbances imparted to the water by a planing 
body give rise to a wave form of motion. on the free sur- 
face, the length of the waves increasing indefinitely with 
increase in the Froude number and being directly propor- 
tional to the latter in the case of the plane or two- 
dimensional problem. Hear the planing surface the general 
picture of the flow as Bhown by tests presents a true jet 
or spray charaoter; i.e., at some distance ahead of the 
body the water surfaoe is practically undisturbed, while 
immediately forward of the body the water is thrown off in 
a spray. 

The ' high-Bpped planing motion of the body gives rise 
to very large accelerations in the fluid and, in this re- 
spect, resembles the phenomenon of impact. The chief 
forces that determine the motion of the particles of fluid 
near the body appear to be the result of the large pres- 
sure gradients. As in the case of impact, it is there- 
fore permissible to neglect the weight of the water. 

The "dynamic reaction of the water is completely de- 
termined by its motion in the immediate neighborhood of 
the planing body. At large Troude numbers the effect of 
the weight BhowB up to any appreciable extent only at 
some distance from the body, so that the flow near the 
body can be considered as part of a flow of an infinitely- 
extending weightless fluid.- The same conclusion can also 
be reached from another point of view. Let us consider a 
series, of motions for which the angles of inclination to 
the' water surface are the same and the wetted portions of 
the bottoms geometrically similar.' Applying the Lagrange 
integral' to the absolute potential motion of the heavy 
fluid,, we obtain the boundary condition of constant pres- 
sure at the free surfaoe In the form 

Report Hoi 252, of- the Central Aero-Hydrodynamical Insti- 
tute, Moscow, 1936. 
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at 2 

where t is the time, " 

<p the' velocity potential, 

v the velooity of the fluid, 

g the acceleration of gravity, 

and y a coordinate taken normal to the initial water 
level. 

Jor steady planing, we have 

at- ax 

whore c is the translational velocity in the forward- 

direction coinciding with the z axis. We Bhall intro- 

.duce nondime'nsional magnitudes with the aid of the rela- 
tions: 

acp^ 3q> x ■ 

- — = o - — ; v = cv x jy ■= hy,, 

where h 1b a certain characteristic dimension. Boundary 
condition (1) then assumes the form 

... ■ . _ |*i + iii + 1| 7l =- 0 • (3 ) 

ox x 2 o 8 
Trom relation (2), it is clear that for the same val- 

C3 

ues of the Troude number T = the motions become dynam 

gh 

loally similar. The values *^L t t., and y. at the 



free surface everywhere except at a very small region at 
the edges of the planing surface are of the order of the 
angle of attaok 0, whioh will he assumed as infinitely 
small in what follows. If the motion of the fluid be de- 
termined by the methods of the theory of waves of small 
amplitude, then in condition (2) it is" necessary to ne- 
glect T i 3 - - At a small value of the Troude number 

7 - the first, and. last terms in (2) are of the order 
gh 
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of 3 and in this case the. weight oannot he- neglected. 

Tor a large value of the 7roudo number, however, I si, 

the third term has the order of 3 a , and therefore to an 
accuracy of the second order of amallneBB. 

|^ « 0 (3) 

Thus, for large values of the Froude number 7, in the 
approximate solution, it is necessary to use boundary con- 
dition (3)., which is equivalent to assuming the fluid as 
weightless* In what follows, the oase will be considered 
where the planing is assumed to take plaoe at high Troude 
numbers. . 

It will bo Bhown below (section it) that in the above 
case the notion of the water at eaoh instant of time may 
be considered as thq result of the simultaneous actio.n' on 
it of a system .of impulses distributed over the. area .swept 
out by thq planing s.urfaoe, the variation with. time of the 
impact pressures onithe free surfaoe being neglected. Thus 
stated, the planing problem differs from that of impact 
on the water only- in that in the former it is necessary to 
take account of the disturbance of the water remaining be- 
hind or, expressed otherwise, it is necessary to takoao- 
count of the asymmetry of the flow in front of and behind 
the body. 

The energy required to maintain the motion of the body 
nay be considered as consisting of the energy of the dis- 
turbed water and the energy dissipated by the dlssipative 
forces of viscosity. The drag may be computed by the 
energy imparted to the fluid by the body in a unit of dis- 
tance. The energy of the disturbance is made up of the 
energy contained in the sprays thrown off ahead of the body 
and that of. . the. 'di.s t urbane es remaining, behind. 

The flow phenomena at the edges in .the case of the . 
two-dimensional problem were made the subj-ect of special- 
study by ifagnar (reference l). In his ,-work, he presents 
methods for the theoretical computation of the portion of 
the drag contributed by . the spray at the forward edge. 
That portion of the. drag due to the disturbance remaining ^ 
behind the body is analogous ' to the induced drag of a wing. 
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Tor establishing the initial disturbance of the main 
flow, the effeot of the viscosity * by comparison with that 
of the inertia of the' water, is negligible. On the other 
hand, the viscosity is of appreciable effeot on the spray 
motion of the water and on the boundary-layer motion. JLt 
small angles of attack, the drag due to the fluid friotion 
at the bottom may be very large. The computation of this 
drag is complicated by the fact that at the forward portion 
the friotional forces on the bottom are directed forward, 
so that the speed of the water in the spray is greater than 
that of the bottom. 

The exact solution of the three-dimensional hydrody- 
namic problem of the planing motion of a body on the Bur- 
face of an incompressible, ideal, and weightless fluid 
presents insurmountable mathematical difficulties and may 
therefore be treated only approximately. The two-dimen- 
sional problem of the steady planing motion belongs to the 
type of flow problems considered by Eirchhoff. The solu- 
tion of the two-dimensional planing problem for a flat 
plate has been given by Chaplygin with the participation 
of Gurevltch and Tanpolsky (reference 2). A comparison of 
this solution with the results of tests shows a very good 
qualitative agreement. Particularly noteworthy is the good 
qualitative agreement of the theoretical law of pressure 
distribution with the pressure distribution determined by 
experiment. The lack of quantitative agreement may be ex- 
plained by the finite span employed in the tests since the 
pressure strongly depends on the span. 

In a fundamental paper on the theory of planing,- 
Wagner, investigating on the one hand flows of the type 
of Eirchhoff and on the other rotational flows about 
thin profiles in an infinite fluid,- showed that for infini- 
tesimal angles of attack, to an aocuraoy of the second 
order of sm'allne'ss, the lift foroe on the planing surface 
is equal to half the lift on a wing of the same profile. 
The flows at the edges of the wing and the planing surface 
are different in character. In particular, in planing, 
thin sprays are obtained at the forward edge, which result 
in an additional drag, whereas, in the case of a wing, 
there are euction forceB -at the corresponding positions. 
According to Wagner, the drag due to the formation of the 
sprays is equal in magnitude to half the suction force on 
the corresponding wing. It is for this reason that the 
drag due to the sprays is not difficult to compute. The 
conclusion also readily follows that the wing always has 
better characteristics than the corresponding planing 



H. A. C. A, .Technical Memorandum- Eo. 942 



6 



-surfaae. Tronr the mathematical point of view, the plan- 
ing .probl-em can easily ^be *e^uc*d^.Ae^we...*haLl-show .below, 
to that of the -wing theory by making use of the approxi- 
mate methodB that have "been applied with great success in 
the theory of thin wings and in the theory of waves of 
small amplitude. The approximations made are equivalent 
to those of Vagner. 

■ 

On figure 1 are given the experimental data derived 
-from -the testa of Sottorf . ( referenoe 3). The tests were 
conducted on flat plates. The aspect ratios* X» l«/b 
where I 1 is the wetted length, as .measured in the tests 
and b the span, laid off on the -x axis. On the ordi- 
nate axis the coefficient \ is laid offx 

Z = (4) 

where -A is the lift foroe 

and p the density of the water, ■ 

Iron these data, corresponding to various conditions 
of motion, it may he seen that the coefficient E depends 

only on the ratio -11. The curve of K = 2 indi- 

b 3(1 + 3X) 

cated by the solid line shows the dependence on the aspeot 
ratio of the value of the coefficient for half the lift 
fordo of the wing as given by the theory of the .lifting 
line for elliptic distribution of the circulation. Compar- 
ing this curve with the tost points, we may observe a 
qualitative agreement. For X - 2, ■ ■ we also have a good 
quantitative agreement. It should be observed that the 
experimental values of the coefficient .K .for the wings 
also do not entirely coincide with. the theoretical values 
and for. small values of X, for example, they differ as 
in the given oase. 

Proceeding to the study of the unsteady or nonuniform 
motion of tho body -within or along the surface of the fluid 

♦In aerodynamics , r the aspect ratio has been defined as b/l. 
When applied to planing* it seemed to us more convenient 
to take as the aspoot ratio the ratio l'/b, since the terms 
I 1 the wetted length and b the- width of step, may be 
considered as definitely established parameters. 
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we are immediately oonfonted with the diff icultiee ohar^ 
aoteristio of these problems. In the case of unsteady 
motion, the mechanioal characteristics of the tody-fluid 
system are not in general defined by very simple geomet- 
ric and kinematic data for the body at the instant con- 
sidered as is the. case for steady motion. If, for exam- 
ple, in a given time interval, we have a constant accel- 
erated forward motion, it would he incorrect to say that 
the hydrodynamio reaotion of the water may he expressed 
as a function of the velocity, acceleration, and geomet- 
ric parameters giving the body position. The reaction of 
the water will, in general', also depend on the character 
of motion of the body in the preceding .time- interval, since 
the body may have previously given a large disturbance to 
the water. 

In the theoretical study of unsteady motion, a num- 
ber of assumptions are usually made for computing the 
hydrodynamio forces. One most often made is the so-called 
"stationary forces" hypothesis. According to the latter, 
it is assumed that during. the unsteady motion the forces 
at each instant coincide with those of a corresponding 
steady motion defined by the same geometric parameters 
giving the position of the body and its trans latioaal and 
angular velocities. The values of the forces determined 
by the stationary forces hypothesis in some instances dif- 
fer considerably from their actual values in the unsteady 
motion. 

Our central problem is to study the relations between 
the forces obtained according to the stationary hypothesis 
and other assumptions and the actual values of these forces. 
Our work constitutes an application to the planing problem 
of the methods and results of the theory of unsteady/ flow 
about a wing. .In addition, we shall attempt to generalize 
and simplify the methods oonsidrerd and present them in a 
form where they may be conveniently used for study and for 
practical applications. 



II. FORMULATION 07 THE PLANING PROBLEM 
AND ITS RELATION TO THE PROBLEM OP A THIN VINO 



We shall consider the main part . of the "f low , where 
the water moveB continuously. The flow will possess an 
infinite velooity at the edges of the planing surface, where 
a thin spray is thrown off. 
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Ve shall determine the ^notion of the fluid on the 
"basis of the following a BMtoptiom i« - - ■- -■ 

'1. Ab shown by experiment ,' there is no -spray sepa- 
rating at the trailing edge at large planing speeds. .The 
water, in this case, /flows off smoothly tangent to the 
"bottom. Ve shall correspondingly require that the velocn 
ity of the fluid at the trailing edge he finite. 

2. Analogous to what is usually assumed in the the- 
ory of a thin wing when formulating the flow conditions 
and in the theory of waves of small amplitude, the "bound- 
ary conditions for the determination of the potential 
flow on the surface of the water and at the bottom, we 
shall transfer along a vertical to the horizontal surface 
coinciding with tho undisturbed water level. Ve shall 
thus reduce- the problem to that of the determination of 
the potential flow of the fluid at the lower portion of 
the half space bounded by the horizontal plane. This ap- 
proximation clearly does not' hold true in the region of 
formation of the spray and at the spray' itself. The mo- 
tion of the water in the spray does not, however, have 
any appreciable effect on the disturbances of the main 
mass of fluid by which disturbances, the chief terms of 
the water reaction, are determined. Vithout too great an 
error, it is permissible to neglect the spray. It may be 
shown that at small angles of attack 0 the thickness of 
the spray and the momentum of the water thrown off in the 
spray in a unit of time are of the order of 0 a (reference 
1). In what follows, we shall assume that 0 is very 
small and shall neglect small quantities of the order of 
0 3 (references 4 and 5). 

3. In the dynamical boundary condition on the free 
surface, we shall neglect the weight of the water and the 
squares of the absolute velocities of the water. It may 
be shown that on the free surface everywhere, with the 
exception of the small region at the edges whore the spray 
is formed, the magnitude of the absolute velocity v is 
of the order of 0 and therefore v 3 ~0 3 . The assumption 
of large planing velocities in- the horizontal direction 
Justifies the neglect of the weight of the water. 

Ve shall sea below that the assumptions enumerated . 
above -immediately reduce the- planing problem to that of 
the -motion of a wing. Tor a steady planing motion, as we 
have seen In the preceding paragraph, such a description 
corresponds to actual physical laws. Tor the unsteady mo- 
tion, this method of formulating the problem when applied 
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to the wing leads to theoretical conclusions that ehov 
good agreement with test results. 

With the aid of the assumptions enumerated above ( 
it is not diffioult to formulate mathematically the 
boundary conditions for the determination of the unsteady 
motion of the fluid. 

Let us consider the two-dimensional problem. In the 
plane of motion the wetted portion of tho bottom profile 
is represented by the curved segment M'H'.ffig. 2). We 
shall study the case where the curve M'ff* differs 
slightlyt in the sense indicated below, from its projec- 
tion MN on the undisturbed surface. 

.Let some point 0 X on the Initial water level be 
the origin of a fixed system of coordinates. The point 
-0, - the center of MIT, we shall take as the origin of 
the movable system of coordinates Oxy, the . y axis be- 
ing direotod vertically upward and the x axis horizon- 
tal. ■ • 

We shall denote by cp(x, y, t) the velocity poten- 
tial of the absolute motion of the fluid, leferring the 
boundary conditions to the x axis, we shall have along 
MN the condition 



3y 



- T, (1) 



where v n is the normal component of the velocity of the 
bottom. She dependence of v n on x and on the time t 
is determined by the geometric and kinematic characteris- 
tics of the planing ' Burfaoe . Let y » f(x) be 'the equa- 
tion of the curve M 1 N ' for- some position of the bottom. 
We shall limit ourselves to the oase where f (x) and 
f'(x) are small. After an infinitesimal displacement of 
the bottom, we shall have for dy/dx, to an accuracy of 
the second order of smallness: 

■ 

^ - P + f '(x) 
dx 

where 0 is the angle of rotation. If o and v are, 
respectively, the horizontal and vertical components of 
the velocity of point M 1 and. u) the angular velocity 
of the bottom, we may write for v n 
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- v -uXa + x) + o 4^ 

dx 



or 

▼n ■ -( T + "») " m z + °0 + °f *( x ) ■ Ti - oux + cf »(x) (2) 
where 

t 4 = -(▼ + u)a). + cp 

Tor a plate forming an angle P with the horizontal, 
we obtain 



(2a) 



If the profile .N'ZT* undergoes bending, then f 1 
is a function not onl7 of . x hut. alBO of the time t. 
In this case, on the right-hand side of formula (2), it 
is necessary to add the normal component of the Telocity 
due to the motion of deformation. 

Sinoe the fluid is at rest at infinity, the Lagrange 
integral on the free surface gives: 



_ + - v 3 + gy = 0 

where g is the acceleration of gravity, and a is a 
coordinate in the fixed system of axes (fig* 2). Discard- 
ing the terms gy and ^ v 3 in accordance with assumption 
3, we obtain 



9cp{a,t) 

at " 



= o 



(4) 



Thus-, on the free' surface <p maintains a constant value 
with respect to time although it may differ in value from 
point -to 'point. We shall assume that, the motion began 
from the state of rest and considering the initial value 
of cp to be zero, we have everywhere on the free surface, 
over which the body has not passed, the condition 



(4a) 



Ve shall denote by S ..the path of the planing body 
on the water surface consisting of the free portion and 
the portion in contact with the bottom at any given in- 
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stant. On .the free portion of the surface S, the func- 
tion cp(a) is determined by the charaoter of the planing 
and may he found from the condition that the velocity of 
the fluid at M is finite. 

The potential of the unsteady motion of the fluid at 
each given instant may he considered as arising from the 
impact of a system of Impulses p^Ca) «= -pcp(a) distribut- 
ed over the surfaces. At a succeeding time instant, this 
surface will he S 1 > S and the system of impulses will 
he p t 1 where, on the.. free portion of the surface S, 

P t ' a Pf 

The potential function w(z) = cp + ]A|r may he extend- 
ed to the upper half piano on the "basis of the Schwarz 
principle of symmetry since to the right of ff on the Ox 
axis, cp ±0.- Let B (fig. 2) he the position of the 
point M when the unsteady motion is set up. We then oh- 
tain the extended function w(s), which is holomorphio 
throughout the plane Bff or AN (the point A denotes 
-<») If the circulation about an infinitely removed con- 
tour is different from zero, and 

cp(x,y) * -cp (x - y) 

>|/(x,y) = >|f(x - y) 

Therefore, at the separating edges above and below, the 
vaiue of cp ' is at eao'h instant the same in magnitude but 
.opposite in sign: 

cp( C) o - q>(Ci) 

It therefore follows that for the motion of an infinite 
fluid determined by the function w(z) the line Bff 1b 
a line of discontinuity of the horizontal component of . 
the velocity of the fluid. The vertical component of the 
veloc ity changes continuously In passing through' Bff. The 
horizontal components above and below are equal in magni- 
tude but oppositely directed. Denoting the discontinuity 
in .the velooity by Y(a), we have 

(d) -« (0) - 1 » (O!) (6) 
da da da 

and . Y(a) may be considered as the density of the vor- 
tices distributed along BM. 
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Let r (a) denote the circulation taken about some 
clotted path.. I|. . (£4g»*. ,2J cutting AJJ at. „ point 0 in the 
direction shown on . the sketch* . ' 

Obvious Xy 

T(a) - -2qX0 4 ) 

whence 

£--*<•> . (•> 

r(<X! ) is the circulation about MH at the instant con- 
sidered. If the circulation about MS is. known as a 
function of time, then, considering the various points on 
AM as the positions of point M at -various instants of 
the motion 



a s» a 0 + / 
^o 



cdt (7) 



and we nay compute Y(a) "by the formula 

dr dt 1 dT 



1{a) = » (8) 

dt da o dt 



Conditions (l), (4), and (4a) and also the extension 
of the flow in the upper half region apply In the same 
form also to the three-dimensional problem with the dif- 
ference only that the discontinuity of the horizontal com- 
ponent of the velocity behind the planing body depends in 
this case not only on the longitudinal coordinate, but 
also on the coordinate in the transverse direction. 

The wetted length M 1 N ' — MIT a 2a in unsteady plan- 
ing may change with time. The wetted length as a funotion 
of the time depends very much on the form and dimensions 
of the bottom, on the amount of immersion, the angle of 
attack, and also on the weight of the water. It is im- 
portant to note that the wetted, length oannot be given as 
a parameter of the motion but is determined in the process 
of solution of the hydrodynamioal problem. 

If a is considered as constant, then, the problem 
we have formulated, coincides ezaotly with that of the 
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determination of the unsteady motion of an Infinite fluid 
about the thin profile M'H 1 vith a surface of disconti- 
nuity separating at the trailing edge* The problem of 
the motion of a flat plate vith a line of discontinuity 
has thus "been formulated by Vagner (reference 6) and that of to? 
oscillations set up has been considered "by Glauert (ref- 
erence 7) and Zeldysh and Lavrentiev (reference 8). Ve 
shall generalise this problem, extending it to the case of 
a variable a and a thin profile of any form. Moreover, 
we shall investigate in greater detail the physical char- 
acter of the hydrodyncmic forces and shall study the forces 
in the unsteady planing motion. 



III. DETERMINATION 07 THE ILUID 7L0V 



To determine the potential function v(z.t) of the 
fluid, we have, on the basis of the assumptions made, the 
following conditions (for notation, see fig. 2): 

1) Outside the segment BH 

4*- = u - iv 

dz 

is a holomorphic and single-valued function 
of z. 

2) At infinity, the fluid is at rest but the circu- 

lation about a contour L a enveloping BH is 
in general different from zero and is equal to 
r o c const. Expressed otherwise, the series 
development of dw/dz near an infinitely re- 
moved point has the following form: 

dw r 0 l c a c 3 

— = + — r + — + . . . 

dz 2tt1 z b b z 3 

3) Along MH on both sides 

dcp 

ay 

and ^ is continuous along BM. 



ay 



4) Along BM, .we have a discontinuity in the hori- 
zontal component of the velocity: 
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The function V(a) Bhould he determined so that 
at point M the velocity of the fluid is 
always finite. 

Let ua consider the function 



T(i) = -A 3 - a 3 .J*. 

as 



Let' us take the hranoh of Vz a - a a , which for 

x > a > 0 giveB a a . Then in approaching MH from 

atoove 

•A 2 - a a -» + i </a a - x a 

and in approaching MU from "below 

A 8 - a a - i ' -y & a - I s " ■ - 

Along AM we have VU a - a 8 ! . 

The function F('z) is holonorphic ..and single-valued 
everywhere outside 'the segment BIT. Sear an infinitely 
removed point the series holow holds' true, namely: 



2-ni s 



Applying the formula of Cauchy to the function r(z), 
we may write 

p(£)dt . i p f(£)dt 

z 

'Li w Ai 



where the integration for hoth integrals is to he performed 
in the counterclockwise direction. 

Contracting the path L x to BH .and expanding path 
L« to infinity, 'we ohtain from;the above expansion for 
F(z) at infinity ■ 

,(.) = JL. pli** * iO) - J(xo - iO) d + 'To. 

2TTi / x Q - z ° 2tt1 
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Denoting "by u x - iv x and u a - iv 3 the values of dw/dz 
in approaching BN from above and below, we have 

and along MH 

v a = v x = -v n 
Vfe tnus find along HH 

T(x 0 + iO) - P(x 0 - 10) = - 2v n (x 0 ) 7a» - x 0 a 
and along BM 

]P(x 0 + 10) - F(x 0 - 10) = 2ui -/x 0 3 - a» = Y(a) Jx 0 * - a a 
Eenco : 

+a _______ 

dw 1 r 2v n -s/a a - x 0 a 

— = m / — dx 0 + 

dz toil J %* - a s J_ a x o ~ * 

a 

i 

1 T Y(q) ^xp 3 - a 3 ^ T 0 
+ — / da + • 

j^i-A 2 - a 2 ^a 0 *o " 2 2TTi/ z 3 - a 2 

(1) 

Between a and x 0l there is the relation 

Xq = cc — a i — a 

The first term of formula (l) gives the velocity field of 
the irrotational motion of the fluid at a given distribu- 
tion of the normal velocities along MN. The second term 
gives the velocity field due to the vortices springing 
from the trailing edge of the plate and distributed along 
BM, in the presence of the fixed plate MN, with density 
V(a). Finally, the third term gives the pure motion of 
circulation about the fixed plate. 

In equation (1) replacing v n by v x - U)x 0 + of '(x 0 ) 
in accordance with formula (2), section II, and performing 
the Integration of the first integral, we obtain 
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■*£ . i Tl (i - - , Jl~. \ ■+ *&■<» - - » V . 

a _ -a 



+a 



ttVz 3 - &a J a - x 0 

-a 

on 

O Y(o)v5c 0 » - a 8 T 0 
+ . -■ / >- da + 



dx. 



2iriVz 8 - a 3 x 0 - a 2nfc/a a - a 8 

o 

whence 



iv^a - -/ 8 8 - a 8 ) - -Jz B - a 8 ) 



8 + 



a p x 0 - i-yA 8 - x 0 8 - z - </a a - aa 

+ £ / f«(x 0 )ln 

" J x 0 + iVka - io a - a - </z 8 - a 3 

-a 



, £i x 0 -A 0 » - a a - a --/a 8 - a 8 

+ — / Y(a)ln —-— da 

^a 0 x 0 +/x 0 b - a 8 - a -Va 8 - a 8 

+ la( , + y B 8 - a 8 ) 8 ( 3 ) 

2tt1 

If the "bottom has the form Of a parabola, wo may put 

y = f(x) *> -ex 8 

whence 

f'(x) = -2ex 
In that case, we shall have 
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A - a \ + l(cu + 2eo) (a - </a a - a 8 ) 8 
^ a/b 8 - a 8 ' 3 -A 3 - a 3 



-A- ■ 1 /* lig&gg - »V ■ r ° (2a 

2111 Va 8 - a 8 / x 0 - a 2nW z a - a a 

^ao 



16 



N.A.C.A. Technical Memorandum No. 942 



(z) = iv x (z - 7z» - a*) - 1(tU + 2ec) (z - </, 3 - a 3 ) a + 



4 



2tt1 / 



Y(a)ln 'o -^o a - a3 -' < „ da + 



a 3 



+ vx 0 a - a a - i - -/a 3 - a 3 



To 



+ — - ln(z + y« a - a 3 ) (3a) 
2tt 4. 

Comparing theBe formulas with formulas (2) and (3), 
ve see that curving the bottom in the form of a parabola 
is analogous to the rotation of a flat plate. (For a 
flat plate we may put f *(x 0 ) =0.) 

The condition that the velocity be finite at point 
M requires that the coefficient of 



2iTi j 



z" - a* 



in the expression for dw/dz become zero when z = -a. 

Thus, +a • a 

^ v n - x 0 3 n 1 y( a ) Jx 0 * - a 3 

- 2 / dx n + / ^ da 0 + T c 



o 



or +a a 



r n - 3 / v n / " ^d i Q . / Y(a) /- ii-^-|da 



Ve shall set 

+1 



a o 



— 1 

so that if 



+ 1 * — : 

i r f «(au) du = 0! (4) 

TT / V 1 + U 



then 



f »(x) 

nss 
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OB 

-5j tea* > iu ak+1 - 'm^ jt*- 8 -^^jiH t n-x 



Pi . . . 

He then have 
+a 



(6) 



Tor a flat plate p x = 0 and v a is the normal ve- 
locity of the plate at a point whose distance is one-fourth 
the vetted length from the trailing edge. 

The condition that the velocity at point M he finite 
In the notation of (5) assumes the form 



• T 0 - 3Trav a » r X Y(a) / ** + * 1 " " da (6) 
Ja Q ^ a, - a 

The above relation Is an Integral equation of the 
Tolterra type of the first kind* Tor given functions 
v a (aj.) and a(a, ) , it is possible with the aid of (6) 
to determine Y(a*. Tor a thin wing, the function v s (a x ) 
completely determines the strength .of the vortices spring- 
ing from the trailing edge in an unsteady motion of trans- 
lation with small angle of attack. 

Tor a steady planing motion with constant velocity c 0 ,. 
we have ti) » V = 0; a « a 0 = const. 

Equation (6) in this case gives: 



T 0 3tfa 0 V ao " 2wa 0 (p 0 + p x ) o 0 

Considering the unsteady motion as a disturbance from 
the state of steady motion, we may write equation (6) thus t 



- anACaVa) --.r'-Vla) / 3a + a x - a ^ (?) 
J »/ a x - a 



a o 
where 

A(av 8 ) i' av 3 - a 0 v 0 ( 
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Tor a rigid wing profile particularly, a and P x 
are independent of time. If, in addition, the horizontal 
velocity at the trailing edge is constant: c ■= const.* 
then 

- 2nA(av a ) = 2na £v+^-c(0-B o )j 

Hence in this case the intensity of the trailing vortices 
does not depend on the ving profile. 

At a » a x the kernel of equation (7) becomes infi- 
nite. To remove this difficulty, ve shall multiply tooth 

sides of (7) "by ■ and integrate with respect to 

ai from a 0 to u. 

We then obtain _ ct x 



t(m) s _ 3n / ,u . r\^r , W /Ji^ ( 

J. J* - a x J J ■/ (a - oOCai - a) 

Changing the order of integration according to the 

Dirichlet formula, we find 

u 



;u) m£ . V(a)] 



. « (u) = / . , Y(a)K(u,a)da (8) 

where 



u ,- 

/ > p / 2a + a x - a 

K(u.a) = / /- r± - da- 

j a v (u - a x ) (a x - a) 

+TT/2 _. 



r 2a + a = g (1 - Bin*) 4* 

"'-TT/a 



(A change in the variables has been made here 

u + a u - a , N 

a, = - sin V ) 

1 2 2 

The kernel K(u,a) is finite and, if a is constant, de- 
pends only on u - a and nay be expressed by elliptic 
functions. The integral oquation (8) may be solved by the 
method of successive approximations. 



da 
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Tor a notion, for, which, the changes in a and- v a nay 
"be" considered av : Infinites inai ( we^obtaiji fr'on .equation (7)* 

Uniting ourBelTOB to values of the first order of anall- 

a s ..... 

ness , 

a. 



- 2n6(av s ) «=• f Y(a) / gft o * a * da 
t/cu V a x - a 



(9) 



Tor a constant, value of . a, equation. (9) coincides 
-with equation (7). Denoting 

'a. 



i. 

Y(a) da = - [r (a x ) - T 0 ] a -ri(tt x ) 

a. 



o 

and setting a 0 + cx x - a = a 0 s, we nay reduce equation (9) 
to the form* 



2TT8(av a )-r«(a 1 )=a 0 J Y( ai +a 0 - »ob)[ y7T"f " *] dB (1D5 



a,-a, 
a o 



If, after a sufficiently large tine interval, the oscilla- 
tions tond to approach a certain stoady oscillating state, 
we nay Unit ourselves. .to...a study of this state. Tor 
steady oscillations we- nay put a 0 = Equation (10) 

then assunes the forn 

2i76 (av a ) - r'taj = a 0 J Y(a x + a 0 - as) £y§±i- 1 j dB (11) 

00 

In equation (10) we cay in genoral consider a 0 = -»»• 
if 6(av s ) is such that for a< a.; 6(av 8 ) = 0. Under 
this condition equation (11) holds true for an arbitrary 
notion. 

IV. THE PROBLEM OT LAUDING OK THE STEP 



Side by side with the problen of stoady planing, it 
is not difficult to obtain a oonplete solution of the 
problem- of the landing with, constant velooity of a flat 
•step' on undisturbed . water; She' solution of this.problen 
iB rendered sinple by the fact that during the innerBion 
with constant velocity in a weightless liquid of a wedge, 
in . the case of the- two-ditoonsional. problen % or .of a cone- 
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shaded body,, iu the case of the three-dimensional problem , 
the motions of the fluid -at different time instants are 
dynamically similar. She problem of the landing on a step 
vlth constant velocity has been considered in the work of 
Vagner (reference l) . We shall consider the problem once 
more here and add a few simple results. 

Let us consider a step which settles vith constant ve- 
locity o on an undisturbed vater surface. Ve Bhall denote 
by 0 . (fig. 3) the angle the bottom makes vith the horizon- 
tal and by k the angle between the horizontal and the 
landing velocity. Ve shall assume the angles 0 and K 
both to be very small. 

Since the motions of the fluid at various instants 
are dynamically similar, the wetted length 2a and the 
circulation along it are proportional to time. Hence 

Y = - — — = const 
c dt 

and therefore equation (7), section III, gives 

- 2Trav 3 - V I / 2a + a i -EH da ' (l) 

J v a x - a 



+ a 

Setting a =» au,— = = s 0 = const > 1,. and replao 

ing v a by c(0 + K), we find 



s o-i 

- 2ttc(3 +k) = y I / B P * 1 " j 

J V B 0 - 1 - U 



du 



o 



whence 



Y 2TT0(P+O ( 2 ) 



^s Q 3 - 1 + ln(s Q + a/s« - 1) 



The minus sign shows that the horizontal velocity be- 
hind the step on the Surface of the water is direoted to- 
ward the .'left. 

• ' 

Vith the aid of formula (2), section III, it is easy 
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to determine .the Vertical component of the volooity of 
: the fluid jpartioieB- oa~the . x« r axis. . Setting,, .x/a_«» b , 
we obtain after a simple computation 



v = o 



2 — (-yr^r/EII 

LJZ^l* ln(B 0 + Jb 0 »-1) X V » - 1 

■ f 

- ln 8B 0 + 1 - J** -l,/s 0 a-lN 1 1 (3) 

B + 8 0 ' .J 

Using this value for the vertical volooity, It is not dif- 
ficult to compute the rise i\ of the water at any point 
with coordinate a* = const 

Ve have : 



a* = cti + a + x 



Si 



nee 



a, i ■ 

ax b ot ; n = — - 



B o " 1 



and 



therefore 



whence 



set 

x = sa 



8 o " 1 



a* = -2 ct 

■o - 1 



a* (so - 1) 

t = 



c(b 0 + s) 



at - - a *<*Q - U ds- 
c(s 0 + s) a 

Thus, for. the rise i\ t .. ve find the formula 
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J 0 0 *4 (•o+-») i 



(4) 



At the forward edge x = a, a » 1, and a* = (s 0 + l)a, 
so that here 



a( 



B ° a I X) f YdB (4a) 



Substituting v from (3) and performing the Integration, 
we find 

n = a(p + k) : , (5) 

ln(a 0 + J* 0 » - 1) + Va o a - 1 

On the other hand, from geometric considerations it follows 
that at the forward edge 

t, = a[3P -(b 0 - 1) k] (6) 

Equating (5) and (6), we obtain the equation for the de- 
termination of b 0 as a function of P/u. 

If I ' = 2a is the wetted length and I is the 
portion of the wetted length below the water level, then 

■Li = , 2ap S 1 ( 7 ) 

I (b 0 - l)ax b 0 - 1 K 11 

On figure 4 is plotted the curve showing how I '/l varies 
with P/k in accordance with formula (?). When p/ — > 00 • 
we have — > oo . It is worth noting that in formula 
(7) c does not enter and therefore the ratio I 1 /! does 
not depend on the landing velocity. 



T. SOLUTION 07 THE INTEGRAL EQUATION TOE STEADY OSCILLATIONS 



*fe shall now consider the solution of the integral 
equation 

3TT5(av 3 ) - r l (a 1 )=a o y^ V(a x + a 0 -a 0 a)^ Jj^r\ " dB <!> 
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for steady oscillations about, a- oertain steady planing 
notion. Est" * ' - ' ■ ■ 

ifct 

6(av a ) = HAe ■ A 0 oos (kt + € ) 

where A = A 0 e ic ; k is the frequency of the oscillations . 

Ye shall conside.r small oscillations, such that the 
amplitude A 0 and the amplitude of the variations in the 
horizontal velocity of the trailing edge are* very small* 
tfe shall, as soune the following law of motion, for the point 
M: 

Cdj_ = C Q.t + Q> J " 

where -a x is a function of time and assumes only small 
values. To an accuracy of the second order of 'smallnesB t 
we obtain ka 

8(av a ) = JLAe °° 

The solution of equation (1) is found by setting 

ka, 
1- 

T»(a 1 ) » EDe 



°o 



whence 

Y( a ) = _-.4E! » £ ["- it "(2) 
da L c 0 J ■ 



Denoting. the abstract number ka 0 /o 0 by n, * wo may write 

ikoi 

ao" * 



Y(a x + a 0 - a 0 s) » - ED ^ e~°o~ e 1 * 4 e" 1 ^ 8 



Substituting the expressions for r (ctj.) and 
V(ai + a 0 - a 0 s) in equation (1), discarding the symbol 

ikoi 

E and dividing by the common factor e °o , we obtain 

2ttA - D = - ijiDe 1 ** J* tT ll **[ J* % \ | - 1 J ds (a) 

The integral on the . righti-hand side, may be expressed by a 
Bessel function, - 7br this purpose, wo shall make use of 
the integral representation of the Hankel functions in the 
Foisson form (reference 9). 
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where 



X = H + icr; n > 0; Xo - f« X x - ^ 



iTor the path of integration, ve may take any curve situ- 
ated in the first quadrant of the v plane. Replacing 

v "by . Ee 1 ^ and assuming that a > 0, we have 

lim /^ e -fi(ffcoB*+^8in*) + iR(nco8*-asin*) (1 _ a lwd ^ B Q 

M "Jo 



and therefore 



/ e (a - 1) 



h£ 1) (M - X p e % ~ ' / a*-( 8 * - 15 " ds (4) 



Using the above formula, we find 



/•*" - i]t. - i^'cx) ♦ £ (a) 

03 

The above relation has been derived on the assumption that 
cr >0 hut it is obviously true also for a ■ 0. Substi- 
tuting 0" «= 0 and replacing i by -i, we find 

*■ oo 

' Substituting the value of the above integral in equation 
(3) , we obtain 

2ttA - - H. i^e^Ciflo 3 ^^) + H^NiOjD 



2 

whence 

46- 1 * 1 A 
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For largo . t&I^ of - rerpL&olng the Eankel functions by 
their asymptotic-expressions Vwe- find - 



Hence 

lim £ - 0 
(i — >• oo A 

Tor small value a of - ix, ve obtain 



4 Ae" 1 ** 



.. . jl + » + £ih ln » 

ti tt Yn 

where 

Y '= 1.781072 

whence 

llm' D = 2ttA 

'. • . ■ - |J — *9 

This value of D may be found if tho circulation 
about MU 1b dotornined in the same manner as in the 
case of steady planing and in the wing theory as in tho 
theory of motion with constant circulation. 

Since the integral equation is linear, it is possible 
with the aid of a Tourier sorles to obtain the solution 
for any periodic law of oscillation. If 

CO 

E ikna 
A n e c o 

then co n= 1 ikna ' 

4 V" 1 A n e- ln ^o— 

** ^LjJ »CHi a, <Bn) - iHi»>(nn)3 

By an analogou"s method with the aid of a Fourier in- 
tegral may be found the solution of integral equation (l) 
at an arbitrary value of 8(av a ), provided 



f*. |8(av a )| da 



remains finite in the region of a considered. 
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71. HYDEOD YNAM I C EOBCES 1ST UNSTEADY PLANIN& 
OB UNSTEADY WING MOTION 



tfe shall now establish general formulas for the hydro- 
dynamic forces vith the aid of the methods given above for 
the determination of the flov. For simplicity, we shall 
limit our study to the case where the bottom is of para- 
bolic form. 

Let Z and Y denoto respectively the horizontal 
and vertical component and H the moment about the center 
of the wetted length of the hydrodynamic force acting on a 
unit width of the planing surface. The pressure on the 
bottom we shall determine with the aid of the Lagrange in- 
tegral. Neglecting the squares of the value of the abso- 
lute velocity of the fluid, we may write 

+a +a 

*P<»»»> dx (1) 
St 



t -y* (p - p 0 h* - - p J 

- - f\j> - Po) fZdx = - P y + p r\ -wMaiil dx 

*J -a J-a. 

+a +a ' 

M -£ (P - Po )xdx = -p^ *^LI*i dx (55 ) 



In the case under consideration, f'(x) = -2ex, so 

that 

I = -P Y + 2e M (I) 

It remains to compute Y and M. The variables a 
and x are connected by the relation. 



a ■ x + f ( c + if) dt + 



const 



and therefore 



d<P(a.t) Sq>(x,t) Sq?(x,t) / . da\ 

at ■ - Tt si - C° + dt; (4) 
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Voting that 



<? (+a) - (?) (-a) = — b -^-g 



where a 



r« = - Y(a)da 



" J ■ 

we obtain o 



a. +a 



-|p(°. + If) r o"|p(» ♦ If) / "M* - p£ <*> 

The value of dcp^ 1 * *j na y be easily found, using 

3t 

formula (3a), Bection III: 

at - at dt d t 2 

a i ' / 

+ 2tt V C+ dt ) (x _ ) 3 - a* « STra-Za^ 



da I" avi _ m^a i + 1 n' t {a) (xi 0 + a 3 ) da "j 

dt L -/a 3 - i 3 2v'a 3 -x 3 + Vaa - x a Vx 3 - a 3 J 

a o 



where u>i denotes <u + 2eo. 

Substituting now ^ from (6) in formula (5) and not- 

at 

ing that 

+a +a 



I ^ — * Lath 

va -■ x 



a 



dx = - tt(x 0 + */x 0 8 - a 3 ) ; 



o -»*o 

a 1 - X o 
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/x Ja.* - x» dx » / ; = 0 

/ - x» 

—a —a 

we find 

T-|p(c + ||)r e + £^^ + p™ Tl ft 



+ pfo + M^i r «<>*<<»>*» P« to £ l v(a) da " 

a 0 - a 0 

To simplify this formula, ve shall make use of integral 
equation (6), section III* from which, noting that 

| x 0 | > a, x 0 < 0 

we have 



a o -o 



*y(a) da 
^x 0 3 - a 3 
a, 



P "Y(a) da 
= 2ira v a - T 0 + a / 

... J ao </x 0 3 .- a a 

Using the above equation, we obtain for the lift foroe the 
final formula 



dt 

+ %- + 2 



* Z' - daN Y < a > da 



Tor steady planing 

V = 0 ; a B const; v x = v a = c0 = const 
and in this case, formula (II) gives 
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T 0 ■ pnao v 8 » pnae 8 3 • 

This Is the well-knoirn formula of Vagner for the lift 
force in planing. 

In studying the motion of the fluid during the im- 
mersion of a slightly inclined 7 "bottom. Wagner identifies 
each of its instantaneous conditions of motion with that 
due to the impaot of a plate of width 2a which, after im- 
pact, has received a velocity equal to that of the sinking 
bottom at the instant of time considered. He thus obtains 
for the lift the formula 



This assumption corresponds in our considerations to the 
substitution in formula (7) of the values c ■ 0; V = 0; 
r Q = 0; aftor which it agrees with the Vagner formula. 
Ve may note that this is the first term on the right-hand 
side of formula (II). This, term takes account of the ef- 
fect of tho so-called added mass. The term 

piTa^c + ^) v a gives the lift force, which may be ob- 
tained on the "stationary hypothesis* 1 and finally the last 
term takes into account the effect of the horizontal veloc- 
ities of the fluid on the free surface, which velocities 
are absent in steady planing. 

Ve shall now compute the magnitude of the moment. On 
the "basis of (3) and (4), we have 

[^n( ot g).,a^ai] te (8 ) 

Jrom formula (2a), seotion III, we find 

x §2. = T , * a + «iL ( a 3 - 2x a ) x + r o x 

bx Va 3 - x 3 2 ./a 3 * x 8 Wa a - x a ' : 

2tt / (x 0 - zV a a - x a 



* Independently of Vagner and on the basis of other consider- 
ations, an analogous formula was established by H. A. Sokolov 
(reference 10) for the dynamic component of the lift in the 
case of steady planing. 
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Substituting in (a) the value of, x|~ from (9) and 

the value of ftP(*i from (7), and noting that 

3t 

+a 



2 



we obtain 

.4 Jin ~ "^i** 



p na* dcu, p -rra da p tibl^ / da\ pa da _ 

M = - r *■ - i — U), — + r Tl ( c + — ) + r r o 

16 dt 4 1 dt .2 x \ dt/ 4 dt 0 

. £S £ /\< a) /fHI da + £2! ( c + flj) / l ^(rt*- 

o °0 

(10) 

Since 



/ V(a) /— — da = T 0 - 2nav a = r 0 - 2irav 1 - Tra 8 u> x 
J *o X ° + & 

formula (10) is considerably simplified to 

prra 4 duu> pna 3 / da\ pa 3 / „da\ / "Y(a)da /_-..,» 

M = - E r v, (c + 2 — )+ r — (c+2 — ) / - (III) 

16 dt 3 \ dt/ 4 V dt/^ o ( y 3C(j a _ ft a 

* . 

Tor a flat plate, u) x is simply, the angular rotational 
velocity w and for a plate bent into the form of an arc 
of a parabola y = -ex 3 , 

U) x = U) + 2eo 

Formulas (II). and (III)' were derived for the case of 
planing. By simply doubling the right-hand sides, we obtain 
the formulas for the lift and moment of a thin wing with 
variable chord in unsteady or non-uniform motion with sepa- 
ration of a line of discontinuity of velocity. 

The values of cp and -f^ above and below are, in 

ot 

fact, equal in magnitude but opposite in sign. Hence, the 
increase in pressure below and the decrease in pressure 

above, due to are exactly equal. Tor this reason, to 

ot ■ 
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obtain the force acting on the wing, it is necessary to 
doubie i the "forces obtained 'f or "the planing motion". .We- 
may note further that for the oase of a wing the term con- 
taining v a in the' formula for the pressure is not neoes- 
sarlly neglected, since the value of the absolute velocity 
of the fluid above and below are ezaotly equal to eaoh 
other, so that this term does not represent a pressure 
Jump in passing through MIT. The term with v* in the 
formula for the pressure on a wing is equivalent to the 
appearance of suction forces directed to the right at the 
sharp leading edge. This concentrated force modifies the 
hydrodynamio action -of the. fluid . as a result of the rare- 
faction which actually take's place at .the rounded leading 
edge of a wing. To obtain the horizontal Component of the 
hydrodynamio force acting on the wing, it is necessary to 
double the horizontal component that acts in planing and 
add to it the suction force. 

tfe shall denote by F the suction force at the lead- 
ing edge of a wing. The value of F is readily computed 
since 

' /(ft)'- 

The integration is performed in the counterclockwise 
direction about an infinitely small circle C with center 

at point z = a. At this point - (*[^) 3 nas tt simple 
pole. Excluding the polo z = a, we obtain 

P = -ott lim (z - a) C^i } 
z — > a n a z / 

dw 

Using the expression'- for — from (2), section III, we find 

dz 

_ pn I" u>a3 o + r f '(x 0 ) ^a a - x o a 

F = t~ v,a - — — - + - - dx 0 

2a L 2 it 7 a - x 0 

+ — r /: da + -s- 

2n / x„ - a 2nJ 

On the basis of equation (6), section III, in the notation 
of (5), (4), section III, we have 
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Therof ore : 

+a a i 

_ . r o r f (x 0 )d» 0 1 n Y(a) da -j 3 , . 

F - 2 ^h + V s 7^=77?^ J^WtT'I (IT) 

particularly for a motion with constant circulation 

t, o r ° + r f l < 3c o) dx o i a 

p - apno h + "X yr^r^J 

tfe Bhall bring out more in detail the significance of each 
of the components of tho lift determined "by formula (II). 

Tor simplicity, we shall consider tho plate as flat; 
Ve then have for tho unsteady planing 
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feast T ,') 

d* V. 2 l J 



▼a. 



T a - pna (o + || ) 

P a / da\ /» x V(a) da 



(T) 



dm 
dt 



yna. 
2 



(o + » w) . 



and for the unsteady motion of a wing with a *» const 



T 2 ■= pirn* 



dv, 

dt 



2pTTQC Yg 

pac 



CCj. _ 

Y(a) da 



Hi ■ ■ 



3TTa 

8 



dUJ o 

IT + pna Vl ° 



(Va) 



P = 2p W a[v 1 + i [ 



■Y(a) da 



da 1 
^5 J 



The forces T x • T Sl and T 3 are perpendicular to the 

plate (fig. 5) , Ti and Y a "being applied at the oenter 

and T 3 at a distance fron the leading edge equal to one- 
fourth the wetted length. 

If we consider the classical problem of the motion 
of a plate of variable width 2a within an infinite fluid 
moving as a potential and continuous flow everywhere except 
at the edge of the plate,, then the general hydrodynamic re- 
action is reduced to suction forces at the edges, the force 
Ti and the moment 



M ' i = 2 M 



) 
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which, we obtain by setting in formula (10) r o = "V = 0 

and doubling the right-hand aide* The force Y x and the 

moment appear ae a result of the so-called addi- 
tional virtual masses. Tor a constant width 2a 

M 7 ! o M x » 2M X 

Thus, the force 1 1 ' and the moment H\ do not depend 
on the circulation T 0 and on the vortices springing from 
the trailing edge. 

If we consider a notion with constant- circulation for 
which motion the condition of flniteness of the velocity 
at point M is satisfied, we obtain the oomponent T x + T a 
and the moment Jt x (reference 11 j. 

The force T 3 or T 3 depends on the strength and 
positions of the separating vortices. As the steady con- 
dition is approached 

*x-> 0, Y 3 -> 0, Y Y 8 

The expression for Y a for unsteady motion coincides with 
the expression for the entire oomponent Y for steady mo- 
tion with constant translational and angular velocities. 
We have, in fact, 

Y a = 2rTpacv 3 (11) 
and, on the other hand, for steady motion, if 

v Q = const; uj = const 

we have (see work of L. I. Sedov (roference 11 ) pp. 30 and 
31) 

y = per 

and 

r = 2nav 8 

that is: 

Y m 2Tfpaov a (12) 
The right-head sides of formulas (ll) and (12) agree 
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not only in outward form but. in actual significance as 
well, if by c and- v a .• In formula Xl2) we understand 
the variable velocities for the given unsteady motion. 
Thus, ? 3 is that value which we should obtain for the 
entire component T if we were to oompute it for the un- 
steady motion according to the rules which hold for 
steady motion, i.e., if we assumed the stationary hypoth- 
esis. 

lor "accelerated motion, if 

a(av a ) > 0 
then V(a) < 0 

i.e., the water breaks away behind at the step. In this 
case, for Js^f-j < c 

always. Hence, the lift force obtained on the stationary 
hypothesis with the virtual additional mass effeot taken 
into account has a value too large. In the case of re- 
tarded motion, the reverse holds, true. Here we are con- 
fronted with on inertia phenomenon for the lift forces. 

The physical explanation given above for the forces 
in formula (V) may bo useful in taking account of the fi- 
nitenoBS of the span of a planing body. Thus, for example 
the component Yi for bodies of various shapes may be ob- 
tained theoretically as well as experimentally. 



VII. FORCES ACTING IN LAND INC ON THE STEP 



We shall determine the hydrodynomic forces acting on 
the step as it is immersed in the liquid with constant 
velocity. * In the notation of section IT, we have: 

tu=0; v x = v a = e(p + k ) 

and 

y 3ttc(3+ k) 

■/b 0 8 - 1 + ln(s 0 + Vs 0 a - 1) 

Formulas (V) now become 
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t ^ : ^dI•r5-. oC^+ •,' t ^ 

I,; = pna (o .+ |£) 0<S 4- k ) 



3 

Bearing in mind tha't ' 

a i + a . ct da c 

— = a - ; a, = ct; . a = r-; — = - — r 

a 1 Bo " 1 dt B o " 1 

and ...... 

• a. 



da. 



— a) 8 - a 3 

we . obtain 



= ■ ln(s 0 + A/a 0 a - 1) 



T x . = P TTC 8 h f _ a 



(»o - D : 



T 3 o. pirc a h 



Ia " pno8h urn** -• 

(- 1 + - k ) i" (b 0 + Dln(B 0 + V* 0 a - 1) 



(s 0 - l)S l. a o a _ i. + i n ( ao + ,/s 0 8 - i) 



where h is the immersion of th-e trailing edge, i.e., 

h a k ct 

Combining equations (1), ; C3) , and (3), we find 

T = piTc 8 he 
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where 

(i + 4) < B o + i)Vi' 0 " - i " " 

£ ^ tL (4) 

<e 0 - 1) 8 CVb 0 » - 1 + ln(s 0 + Vb 0 » - 1)3 

Tfo have further 

Mi - f <*i + *»> 



Thus, in landing on the atep with constant velocity the 
general hydrodynamio reaction reduces itself to a single 
force applied at a distance from the leading edge equal 
to one-fourth the wetted length (fig. 6). 

From formula (4) , it is clear that this force de- 
pends essentially on the ratio 0/ k. Figure 7 shows the 
naturo of this dependence graphically. It is interesting 
to observe that the curve for C has a minimum, a faot 
whioh indicates tho existence of on optimum ratio 0/ k 
in landing. 

The impact force computed from formula (4) will ac- 
tually he larger than the force applied at the step since 
we considered a flat bottom of infinite span. The finito- 
ness of the span and the presonoe of a 7 angle result in 
a decreased water reaction. 



In oonneotion with the physical explanation of the 
hydrodynamio forces given in the preceding paragraph, it 
is interesting to compare with one another the values of 
I 1( T 8 , and T 3 for the case under consideration. Since 
all these forces are proportional to time, their ratios 
will he independent of time. For comparison, their varia- 
tion with tho parameter 0/ k -are shown graphically on 
figure 8, the values Yj/k, T a /k, and T, /k , whore 

k a P™(P ♦ n)o 8 t 
(so - D a 

being laid off on tho ordinate axis. 
Trom formulas (1), (3), and (3), 
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(b 0 + 1) ln(B 0 + ^b 0 « - 1) 
k -/b 0 » - 1 + ln(s 0 + Jb 0 * - 1) 

We see that for small values of £ ( T 3 la comparable 

with T a and in this case it is therefore inoorreet to 
neglect to take into account the effect of the tangential 
velocity of the vater bohind the stop. 

The limiting eases as — — > » and — — >0 are 

K K 

shown on figures 9 and ,10, respectively. In the first 
case, vo have 

lln — * ■ lim 8 0 - 15.645; lim — a - -2.995 
k k 

in the second I 

lim la ■ l| 11a Zp. s -l 

*t XL 

Till. TOBCES ACTIHG DUBIJTG THE STEADY OSCILLATIOH 01 A WING* 



Wo shall consider the magnitude of the hydrodynamic 
forces during the steady harmonic oscillation of the ving. 
Wherever magnitudes appoar in complex form, it is to he 
understood that only their real parts are considered* 

Let the oscillations he of tho form 

5(av a ) ■ Ae lkt and 8(av 1 ) ■ Be 1 ** 

where 

A = A 0 e 1(6 i + C «>| B = B 0 e i€ * 

A 0 and B Q boing Bmall values. In the case of the wing, 
a s const. 

Considering, for simplicity, o to he constant, we 
find from formulas (Ta) 



•The prohlem of the steady oscillations of a plate has heon 
treated by other methods In the paper by Glauert (reference 
7) and by Eeldysh and Lavrentiev (reference 8). 
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?i «.prraikBo ikt ■■ (D 

*a - T sp - 2pTT0Ae lkt (3) 
From formulas (2) and (6), section V, we have 

Y(a) = — — (3) 

• CHj 8 >(|*) - iH< 8 >(n)] 

Using the above and setting 

x 0 = a - a, - a = -as 

we find 



*3 



°1 - lk0L l ■ x * _ 4 fl 

/» 7(a) da _ 4pciAe'Q P e" da 



The integral on the right-hand side is expressed by Hankel 
functions of the second kind. Prom formulas (4), section 
V, at p = 0, and replaoing i by -i, we obtain 



1 , V 

e ds Tri (a) 

— Ho ^ 



*/a a - 1 3 

oo 

Ve shall furthermore denote 



Hp 8) (n) , , i X (n) 

-r-r r-r *v (lOe AXM '-. (4) 

H 0 3) („) - i Hl (a) (n). 
Ve then have finally 

T 3 - -3pTTcA V e l(l,:t+X) (5) 

The magnitude V as sben from formulas (3) and (5) iB the 
ratio of the amplitude of the lift component, arising from 
the presence of the . trailing vortices , to the increase in 
the lift force computed. in accordance to the stationary 
hypothesis! .. 

max T s 
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ka 

On figure 11 is. shown how v depends on \x = — . When 

1 ° 
\i «= 0, v a 0 and when 11 = 00, v ■= — . 

2 

The angle x is the phase shift of the foroe - T 3 
with respect to the foroe T a - T ao . Figure 12 shows the 
depondenoe of X on V" When \i = 0 

and when n = », 

X = 0 

The moment M x in the case considered is now oppressed 
as follows: 

M x - M 10 = - — a (A - B)il=e ikt + pir acBe ikt (6) 

4 

ka 

At large frequency, i.o., when — > 1, the effect 

c 

of the trailing vortices results in a decrease in the lift 
force to one-half its value obtained on the stationary 
hypothesis. 

* 3 = - I <T B - T ao ) 
ka 

At small valuos of >0, i.o., at large horizon- 

c 

tal velocities, the value of v is near zero and there- 
fore Y 3 — 5*0. In the latter case, tho force obtained 
according to the theory of motion with constant circulation 
is near in value to the general actual force. It should 
be observed that, as may bo seen from figure 12, Y 3 ap- 
proaches zero rather slowly when p, — > 0. 

The suction force P is determined with the aid of 
formulas (7a). Since 

_1_ p 1 V(a)da = Y 3 5-s a vc i(kt+X) 
2n / </s 0 3 - a 3 2-rrp ac ~ a 

therefore 

1 3 



_ Bon f ^ . i(kt+X)T 

P = -J- |^av 10 + 8(av 1 ) - EAU e 'J 



whence 
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P «= ^£H + B 0 cos(kt + e x ) - A 0 vcos(kt + € x + e a +X)J 



t t 



In conclusion,' we shall consider the projection of 
the.- general hydrodynamlo foroe on .the horizontal direc- 
tion. .Vhen the ving has the form of a flat plate t ve have 

I - -(Tx + Y 8 : 4 + T 3 )0 :+ P 

where. 0 Is the small angld of Inclination of the plate 
to the horizontal I 

= P 0 + j^j A o ain(kt + e x + e a ) - ein(kt + 6^) 

Tho mean value of the horizontal component over a. 
! period Is . 

STT/k 

*° = 2?F /" T dt B Tt{ a ° Bi, ( v " 2 a¥ 8ln *) + AoBoToob e a 

- 2v cos <X + e a ) + 2 ^(wsln (X + e a ) - sin c 3 )j}< 8 > 
Tor translatlonal oscillations 
Ao " B o'«. € a '= 0 

In this case : 

P" -o . , . ,.a _ „ coa x |> o 



£i + v 8 - 2 u cos xj: 



a 

l.e\, we have' a forward thrust on the wing.. 

Tor purely rotational oscillations about the center 
-of the plato, B 0 = 0 and therefore 

. ' jc n .i t ■ P" Aq 8 , v [i> - .2 JL sin X j 
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IX. FOECES DUH I HG STEADY QSCILLATJOH OF A PLAHIHG FLATS 



In applying the general formulas derived in eeotion 
VII, it ia first necessary to detern}.no- a and 

4r Where a = —• is half the wetted length. As tests 

at <o 

have shown, the wetted length I 1 is always somewhat 
larger than the part I . under the surface level. On 
figure 13 are presented data obtained from tests of 
Sottorf (reference 3) and those conducted at the Central 
Aero-Hydr odynamical Institute by Perelmuter (reference 
12). The ratios are plotted as ordlnates against 

the aspect ratios Jl <= A as abscissas. The tests were 

b 

carried out on flat plates. It is readily observed from 
consideration of the plotted data that 1* > I. At large 
aspect ratios 

I ' - 1.1 \ 

At small aspect ratios corresponding to a largo span I* 
may exceed 1 considerably. 

To determine theoretically the ratio for the 

case of steady planing motion 'of a plate of Infinite span 
is not possible if we limit ourselves to the consideration 
of the motion of a weightless ideal fluid. Actually in 
the flow about an inclined plate the height h of the 
trailing edge above the level of the freo surface (fig. 
14) approaches infinity as the distance from the plate in- 
creases (reference l). For this reason, it is also impos- 
sible from such considerations to determine the wetted 
length below the surface level. In the flow of a heavy 
fluid with spray formation, the fluid in front of the 
plate cannot rise higher than a certain level, depending 
on the planing speed. From the test data presented,, it 
is seen that the wettod length depends essentially on the 
span of the planing plate. 

Let us make the assumption that, in unsteady motion, 
\ ' depends only on I and that 

I • = nl (1) 

where n is a c ons'tant ' coefficient . The value of n may 
be obtained from test data referring to that stationary 
state with respect to which we consider the disturbed mo- 
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tion to take place. Glauert and Perring in their paper, on 
the stability of . a_ seaplane in planing assume l 1 = I, 
i.e., n ■ 1 (reference 13). 

It is necessary here to make the following reserva- 
tion. On figure 13, it is seen that for small aspect 
ratios n- dopends very much on the aspect ratio. With 
unsteady planing,- the aspect ratio is variable, so that 
for those cases where it is small we have no previous as- 
surance that it is always permissible to neglect the vari- 
ations' in n; At large aspect ratios, however, the varia*- 
tions in n are negligible. 

Considering n as constant and the disturbances 
small, 8a may easily be obtained, tfo have 

25a = 61 » = n8l = n8 - = - ^fil 

0 Po Po 

where h is the immersion of the trailing edge, y «= 8h, 
6 = 60. 

••-Wr# ■ <*> 

— _ n dy _ a 0 d.9 . . 

dt " 20o at 3 0 dt 

Ve shall now proceed to the determination of the hy- 
drodynamic forces. Retaining in formula II only small 
quantities of the first order, we have 

Pn a 0 d(avi) 
™ = ~ dt~~ 



+ pTT.o 0 8(av a ) + pn a o c o 0 o (8 o + |> 




•Y(a>aa _ ( 7I ) 
Jx 0 * - as 

If M a 1b the moment of the hydrodynamic forces with 
respoct to the trailing edge, then 

M a = M + (a 0 + 8 a)(T - Y 0 ) + (a 0 + 8 a)Y 0 

From (III) and (VI), to an accuracy of the second order 
of smallness, we find 
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M a - M a p = - ELSfil i£ + BH^l iSSil + ££f>°° 8 ( aTl ) + pnaccos (av 8 ) 
■ 11 . lo at 3- at -2 

3a 

+fma 0 Co 3 P 0 (|8c + | 5ft) + | pTra 0 o 0 a (3 0 6a + -g 2 " T 5 (Til) 

The variations 6 ) and 8{av 3 ) nay he expressed fcy 
the klnenatlc elements of motion of the trailing edge. We 
have 

. .a 

av 3 = (a 0 + 6a) (c 0 + 6c) (P 0 + 8) + (a 0 + 8a) |£- a °* ^ ^ 

and fron formula (?), retaining only small quantities of 
the first order, we ohtain 

8(av 3 ) = u + a 0 || - ^ ^ + a oPo 6c (4) 



and similarly 

sCavJ = *f y + a 0 || - a 0 a || + a 0 P 0 8o " ■ (5) 



In determining the forces, the greatest difficulty 
met with is in the computation of I 3 . To obtain this compo- 
mont, we shall consider small steady oscillations for which 

8(av 3 ) = H Ao ikt (6) 

The integral equation for the determination of V(a) in 
this case is the same as for the wing and therefore 

*3 = 1*3 = " HPTTACVO 1 ^*^) (7) 

Ve shall now consider in detail the forces acting on 
the plate moving according to the law 



y = E Pe ikt 
6 = fi Qe 
8c = H ik ffo ikt 

where F,- Q, and H are -complex constants. 



(8.) 
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Using WBlatrone'i '(3), (3),. (4)., (5), and (-8).. wo find 
f r on f omul as . . ( 7.1 ), and _( V I i ) : 

Y - Y 0 = He^ (Elil n(l - Vf} 1 *)? - prrao(i + n - i?e iX )ikP 



- Z0*L k» P - fm.»o(a - f e**)lkq + k a q 

i 

+ pTioaP(2 - ue 1X )ikN - P.k 3 S } (9) 

and 

M B - M ao - ^"{f P"^n(l - |ve iX )P 

+ | prra a c(l + n - ve iX )ikP - £^?-kaF . £ pTTa a c a <4 



- | pttb=c (2 I - -| Ve lX )ik(i + ^ prra*k a Q 

+ I prra a oP(2 - ve 1X )ikN - f~.pk a n} ClO) 



In the above formulas, the subscript o has "been omitted 
from those magnitudos which correspond to steady planing. 

Toflrmulas . (9)- and (10) wore- derived "by us from the 
general expressions for the hydrodynamic forces in un- 
steady planing. It is interesting to compare them with 
the results which may he obtained on the stationary hypoth- 
esis or on the assumptions of Glauert and Ferring. The 
stationary hypothosis leads to the formulas following: 

T = T a = prra^o + v a 

where M x is' the moment with respect to the center of the 
wotted length. The moment with respect to the trailing 
edge ( is equal to 

. ■ 

M 3 = M, + aY 3 = -EUS 3 . ( c + 3 J^W + pTTa 3 Co + Ta 
' 2 ■ ^ dt/ 1 > dt / 
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Batalning snail ,te.rns of. tji'e first ardor only 

T * T Q = p-Te6(av a .) + pTTaoP^6c..+ ^ 
H a - Hjo b i2E£2. 8(av,) + p-rac6(av a ) + -J p-Tac a p8a 

+ £ pira a cp ffic + 1 ^ 
2 K \ 3 dt/ 

The above fornulas applied to tho harnonio oscillations 
under consideration give 

i 

T. - T 0 - 1 o lkt {£2|fE'p + prrac (l + |) ikP 

- § piTa a cikQ + 2p-raopikN J- (11) 
and- ■ 

Ms - M ao = j£e lkt pnac. a nP + £ pna a c (l + ■^■^ ikP 

- | pira a c a q - 3 p-Ta 3 cikQ + 3p-Ta a cpikN}' (12) 

• ■ , 

Olauort and Perring, in thoir papor cited above, 
nake the following assumptions to obtain tho hydrodynanic 
forces x 

■ * 

: l) In 1 unsteady planing notion, the .Irif.t force is con* 
• putod according to the foruula 

T m qev a p 

where q is a constant coefficient proportional to the 

width of the plate and to the density of the 
fluid, 

s the wotted lon^th below the undisturbed water 
level (fig. 15) , 

v the magnitude of tho velocity of the point of 

the plate coinciding with tho center of pros- 
sure, . .' • v. 

P tho angle between the plate and the velocity 
direction. 
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2) In unsteady planing riotion,' the force T is ap- 
plied at a distance cp B from the trailing edge 

whore ~<p is the value , in steady planing. Tor 

a plate <p » 3/4. 

Using these assumptions f discarding terns of higher , 
order of smallness, and setting s = 2a; v 0 ? c, vo ob- 
tain - 

Y - Y «= _ 

° 2dP 




M a - M ao = ^ y + — ~ 

The values of the hydrodynamic forces obtained for the 
harnonic oscillations of a plate on the basis of the vari- 
ous, assunptions made are compared in the fable belov. -In 
transferring the data fron formulas (9), (10), (11), and 
(12) to the table, the magnitude pTrac a 0 has been replaced 
by T 0 . 



Translation by S. 
national Advisory 
for Aeronaut Icb. 
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Committee 
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Comparison of Hydrodynamlc Forces Obtained on the Basis of 
Various Assumptions for Steady Harmonic Oscillations 
with Respect to the Stationary Planing Condition 





In the formulas for the lift T - X 0 


Coeffi- 
cients . 


General 
formula 

U.-*S 

c 


Slow 
oscilla- 
tions, 
large 
horizon- 
tal vel. 

v=o 


Rapid 
oscil- 
lations 

y-i- 


General 
formula 
accord- 
ing to 
"sta- 
tionary 
"hypoth- 
eses" 


Accord- 
ing to, 
assump- 
tions of 
Glauert 

and 
Perrihg 






. n 
2p 


,71 

: ■ 4|3 ~: 


n 

r? 


'i 

2p 


y 0 Pe»* 
a 


- P 


(1 + n) ^ 


V T 27 p 


U ^ V 2 


. ¥ 
P 






ft 
2P 


2P 


0 


0 




0 


0 


0 


0 


0 


y« Qe"" 


-(2-|VX)f 


2lf 
P 


3 <|t 
4,3 


2 p 


1 1 '1* 

2 p 




2? 


ft 

2? 


2P ■ 


0 


0 


Y^Ne"" 
a 


(2 — ve ''■) If- 

ft 
2 


2lf 

ft 

. 2 


2 


2i> 
0 


2(> 
0 



In the formulas for the moment with respect to 
the trailing edge: M2 - M20 
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Figure 2 
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Pigs. 3,4,5,5,7 
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Figs. 8,9,10,11,12,15 





Figure 15. 
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Tigs. 13,14 
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